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Abstract. We study percolation on self-dual hypergraphs that contain hyperedges 
with four bounding vertices, or “four-edges”, using three different generators, each 
containing bonds or sites with three distinct probabilities p, r, and t connecting 
the four vertices. We find explicit values of these probabilities that satisfy the self¬ 
duality conditions discussed by Bollobas and Riordan. This demonstrates that explicit 
solutions of the self-duality conditions can be found using generators containing bonds 
and sites with independent probabilities. These solutions also provide new examples 
of lattices where exact percolation critical points are known. One of the generators 
exhibits three distinct criticality solutions (p, r, t). We carry out Monte-Carlo 
simulations of two of the generators on two different hypergraphs to confirm the critical 
values. For the case of the hypergraph and uniform generator studied by Wierman et 
al., we also determine the threshold p = 0.441374 =b 0.000001, which falls within the 
tight bounds that they derived. Furthermore, we consider a generator in which all 
or none of the vertices can connect, and find a soluble inhomogeneous percolation 
system that interpolates between site percolation on the union-jack lattice and bond 
percolation on the square lattice. 
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Percolation is a fundamental model in both statistical physics and mathematics 
[1, (2, [3, H El i;6j. It is concerned with the formation of long-range connectivity which 
occurs when the occupation of sites or bonds exceeds a critical threshold. Finding 
the exact threshold for different lattices is one of the primary goals in the study of 
percolation. Precise thresholds are necessary for applying percolation models to real 
systems, and for studying the behavior of systems near the critical point. For many 
years, exact values of threshold had been known for only a handful of lattices (e.g., 
square, triangular, honeycomb, kagome, and bow-tie lattices, for site and/or bond 
percolation) [7, [8, Recently, percolation thresholds have been found for a broad 
class of lattices that can be represented in the form of three-hypergraphs self-dual 
under the triangle-triangle transformation mm- In a hypergraph representation, 
the hyperedges (in the shape of triangles, squares, etc.) can represent any collection of 
bonds and internal sites, including correlated sites and bonds; all that matters is the 
connection probabilities between the boundary vertices. For a triangular generator, 
there is a single nontrivial self-duality condition which gives a unique percolation 
threshold mm na m 

Prob(all vertices connect) = Prob(none connect) (1) 

and this has led to many exactly soluble lattices when applied on self-dual three- 
hypergraphs D23 However, several important lattices correspond to hypergraphs with 
edges of four vertices (four-edges) and cannot be represented on a three-hypergraph, 
and it is desirable to find thresholds for such lattices as well. Bollobas and Riordan 
m discussed the self-duality conditions for self-dual four-hypergraphs, and recently 
Wierman et al. [15J used those conditions and a stochastic ordering method to find 
bounds for a 16-bond uniform probability generator on a certain four-hypergraph. 

The purpose of this paper is to find explicit self-dual four-vertex generators and 
to confirm with Monte-Carlo simulation that when they are placed in a self-dual 
hypergraph containing four-edges, the system is indeed at the threshold. In contrast to 
the case of hypergraphs with triangular generators, the conditions for a square do not 
reduce to one single equation but to three as we will see below; therefore, we allow for 
three different probabilities to ensure that the set of equations have nontrivial solutions. 
It will turn out that the number of solutions depends on the generator we choose. We 
study the hypergraph introduced by Bollobas and Riordan (hypergraph A) shown in 
Fig.0 along with another hypergraph based upon the (3 2 ,4,3,4) Archimedean lattice, 
which we call hypergraph B, shown in Fig. [2} We consider in detail two generators, 
one with 12 bonds (generator I), and one with 16 bonds introduced by Wierman et al. 
[15] (generator II), both containing three distinct bond probabilities p, r and t as shown 
in Fig. [3| We verify the criticality of these solutions numerically. We also study the 
uniform case p = r = t for generator II on hypergraph A numerically and find that the 
threshold falls well within the bounds found in m Furthermore, a variant of generator 
I (which we call generator III) is considered by replacing the inner square in Fig. [3ji 
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with a single site with probability t as shown in Fig. [3|:. We do not however, perform 
Monte Carlo simulations for this case as this generator is not the primary focus of this 
paper. 



Figure 1 . (a) Hypergraph A, and (b) showing the self-duality. 

Hypergraph A is a four-uniform self-dual hypergraph—i.e., consists of only four- 
edges connecting boundary vertices in a self-dual configuration, as shown in Fig. [lj In 
constructing the dual, a vertex is put in each empty polygon (or face of the hypergraph), 
and a dual hyperedge is drawn around each hyperedge. Hypergraph B is a bit different 
from hypergraph A in that it also contains two-edges (ordinary bonds). These bonds 
have probabilities pi and 1 — pi in a manner that leaves the hypergraph self-dual (Fig. 
ID' Because the hypergraph is self-dual, the value of pi does not matter at the critical 
point, as long as the generators in the four-edges satisfy the duality conditions. We can 
therefore manipulate pi and nothing should change. A variant of this hypergraph which 
only involves four-edges is obtained by setting p\ = 0 (Fig. [2^). 



(a) (b) 


Figure 2. (a) Hypergraph B with four-edges (squares) and single bonds of probability 
Pi and 1 — pi. (b) Hypergraph B with probability p\ set to zero. 

We note that another approach to finding thresholds for some lattices has recently 
been put forth by Grimmett and Manolescu [IS] , for lattices that can be represented in 
an isoradial form in which each polygon can be inscribed in a circle of equal radius. This 
method can be used to find a geometrical proof HZ! of Wu’s criticality condition for the 
square checkerboard lattice [T8j, however, it cannot be used to study the types of square 
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Figure 3. (a) Generator I (b) Generator II (c) Generator III 

lattices considered here. We also note that methods have recently been developed to 
find approximate thresholds for many lattices to extraordinarily high precision [191 f20] . 

In the following sections, we discuss the definitions (Section 2), the general self- 
duality conditions (Section 3), the derivation of the explicit critical points of the two 
generators (Section 4), Monte-Carlo simulations (Section 5), the analysis of the critical 
manifold for generator I (Section 6), the derivation of a generalized union-jack lattice 
for site percolation (Section 7), and conclusions (Section 8). Explicit polynomials for 
the generators are given in the Supplementary Material. 

2. Definitions 

Consider a generator G with four boundary vertices A, B , C, and D. Denote the dual 
generator by G*. 

Any configuration (i.e., designation of internal edges as open or closed) on G determines 


c 



Figure 4. Solid lines represent a four-hyperedge with boundary vertices A, B , ( 7 , and 
D. Dashed lines represent the dual hyperedge, with its boundary vertices A*, £>*, ( 7 *, 
and D*. 

a partition of the boundary vertices into clusters of vertices that are connected by 
edges. We use the following notation defined in m : connected vertices are grouped into 
clusters, and distinct clusters are separated by a vertical bar. For instance, P°[AB\CD] 
means that A and B are connected, C and D are also connected, but the two sets are 
not connected to each other. We also introduce the quantities P n for isotropic systems 
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Pi = none of the vertices connected = P G [A\B\C\D\ 

P 2 = only two nearest vertices connected 

= P g [AB\C\D] or P g [BC\A\D} or P G [CD\A\B] or P G [DA\B\C] 

P 3 = three out of four vertices connected 

= P g [ABC\D] or P g [ABD\C] or P G [ACD\B] or P G [BCD\A] 

P 4 = all four vertices connected = P c [ABCD] 

P 5 = only two vertices connected diagonally 
= P g [AC\B\D] or P g [BD\A\C] 

Pq = two unconnected pairs = P°[AB\CD] or P G [AD\BC] ( 2 ) 

Note that we cannot have P G [AC\BD]\ otherwise the graph will be non-planar. 
Normalization requires 

P\ + 4 P 2 + 4 P 3 + P 4 + 2 P 5 + 2 P 6 = 1 (3) 

3. Self-duality conditions 

By the duality relationship between G and G*, an edge in G* crossing an open edge in 
G cannot be open. Immediately it follows that P G [ABCD] = P G *[A*\B*\C*\D*] (Fig. 
|2J). Assuming that at the critical point, G and G* should have the same probability of 
connection between the vertices, we get the first self-duality condition which needs to 
be satisfied: 

P*(p,r,t) = Pi(p,r,t) (4) 

analogous to Q for the triangular case. However, here we get two additional conditions: 

P 2 (p, r,t) = P z (p,r,t) (5) 

Ps(p, r,t) = P 6 (p,r,t) (6) 

as given in [14J. Figures [5] and [ 6 ] illustrate these two additional duality conditions. 

The three relations yield three nontrivial equations, unlike in the case of hypergraphs 
with triangular generators where there is only one nontrivial equation ([l] for criticality. 
It is for this reason that we chose three distinct bond probabilities (p, r, t) within the 
generators instead of only one uniform probability p. 

4 . Derivation of the critical points 

To find the critical points, we need explicit expressions for the various P n (p,r,t). For 
the two generators, the six probabilities defined in (§ will be of the form 

N 4 4 

P n (j,,s,t) = J2'EYi c%h > i q N -Vs l -H k u i - k (7) 

z=0 j =0 k =0 

where q = 1 — p, s = 1 — r, and u = 1 — t, and N = 4 for generator I and N = 8 for 
generator II. 
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c* 



Figure 5. If A, B, and D are connected in G (a case of F ’>), only A* and B* can be 
connected in G* (a case of P 2 ). 


C* 



A* 

Figure 6. If A and C are connected diagonally in G (a case of Pf), then B* and C*, 
and A* and D* can be connected in G* (a case of Pq). 

To find the c[jl, we used the method of exact enumeration, as follows: Go through 
every possibility of bond placement (there are 2 12 and 2 16 possibilities for generators I 
and II respectively); calculate the probability of percolation for each bond configuration 
noting that each configuration corresponds to a unique set of {i,j, k}] then count all 
the possible configurations that give the same probability of percolation (i.e., the same 
{i,j, k})] also keep track of the vertex connectivity for each configuration using standard 
cluster algorithms. The last two steps identify for a given n and {/', j, k}. The results 
of these calculations are given in the Supplementary Material. 

Having found the probabilities, we then solve and (|6j) numerically to find 

the self-dual points. For generator II we find a single physically meaningful solution: 

p = 0.50003748516983960964249978886157073806 
r = 0.19560520467878219513987817879753912889 (8) 

t = 0.45118455132998743413937451492248773536 
while for generator I, we find three distinct solutions: 
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First solution: 

p = 0.07878784860198622232682114813885158623 
r = 0.72008740215372047101169734956905056064 

t = 0.56558686100571432867919457371281967299 
Second solution: 

p = 0.14732762147606095852100839931621084702 
r = 0.59819060855262599710283661843932567345 

t = 0.64087266343391141007369901653243921773 
Third solution: 

p = 0.18026397627734291307973377345301879693 
r = 0.50904978773535203861705700693678065244 
t = 0.75916858817391565479639455248414312467 
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( 9 ) 


( 10 ) 


( 11 ) 


As we will discuss in Sec. [6| these three points lie on a manifold in (p, r, t ) space. By 
slightly changing the values we can find the other points on the critical manifold, which 
now depends upon the hypergraph being considered. We find for a range of probabilities 
on the critical manifold away from these three exact solutions, the duality relations are 
very closely (but not exactly) satisfied. 

We also look at generator III (a variant of generator I), in which the vertices on 
the inner square are correlated in such a way that either all of them are connected 
with probability t or none are connected with probability 1 — t, i.e. the inner square 
in generator I is squeezed into a site that is open with probability t or closed with 
probability 1 — t (Fig. |3j). In this case we find a single physical solution: 


p = 0.19738202171710149917476592797778709312 
r = 0.44960772662591992436180558214955545825 (12) 

t = 0.99956609784920984836488846372795998038 


Interestingly, the solution for t is very close to 1 but not exactly 1. If t = 1 were a 
solution, then we would have a two-probability generator that satisfies the three duality 
relations. Likewise, a generator with four outside bonds and two crossing diagonal bonds 
does not satisfy the duality conditions. The values of p, r, and t above are closest to 
the third solution to generator i@, but still quite different (especially t ). 

For comparison, we present Pi(p,r,t ) for the solutions given in equations ([8-12 ) in 
Table [l} We observe that the values in any of the columns are very close to each other. 
In fact if we plot for example P, vs. P 3 they nearly fall on a straight line. In Sec. [6] we 
give a plot of Po vs. P\. 




Percolation on hypergraphs with four-edges 


Table 1 . The values of P,(p. r, t) for the self-dual points given in (js]-|T2j). 


generator 

Pi (= p i) 

A (= A) 

A (= A) 

I (soln 1) 

0.164163083 

0.071552907 

0.024812643 

I (soln 2) 

0.164308779 

0.071547560 

0.024750489 

I (soln 3) 

0.163414538 

0.071598550 

0.025095631 

II 

0.166459728 

0.071431975 

0.023906186 

III 

0.162628061 

0.071642238 

0.025401492 


5. Monte-Carlo simulations 

We carried out Monte-Carlo simulations using a Leath-type of growth algorithm on a 
lattice of size 16384 x 16384; clusters were grown up to a size cutoff 2 20 = 1048576, 
guaranteeing that the boundaries of the system were never reached. This allowed us 
to find an unbiased estimate of P> s = the probability that a point belongs to a cluster 
whose size is greater than s. According to scaling theory, we have at p c , P> s ~ s 2 ~ T 
where in 2d, r — 2 = 5/91, and for p close to p c , we have 

s t ~ 2 P> s ~ af(b(p - Pc)s a ) « A + B(p - p c )s a (13) 

where f(z ) is the universal scaling function, and the last term above follows from a 
Taylor series expansion of f(z ), where o, b, A, and B are non-universal constants which 
are specific to the system being studied, while a = 36/91 and f(z ) are universal. Thus a 
plot of s t ~ 2 P> s vs. ( p—p c )s° yields a straight line (for large s ) with a slope proportional 
to p — pc . When p = p c , s T ~ 2 P> s is constant, apart from deviations for small s where 
scaling is not valid. 

5.1. Generator I on hypergraph A 

In this case, because we have two parallel bonds where the square generators touch, 
it was convenient to replace those two bonds having probability p by one bond having 
probability p' = 2p — p 2 [9], as shown in Fig. [7| Figure [8] shows the results for the 
simulation of 250 000 samples on a lattice of 16384 x 16384 for generator I on hypergraph 
A, at the predicted critical point of the first solution Q and with values of p equal 
to 0.0001 above and below the critical value 0.07878785.... It can be seen that Q 
corresponds to a critical point within high numerical accuracy (at least 10 -5 ). Likewise, 
we verified that the other two solutions are also at the critical point within this error. 

5.2. Generator II on hypergraph A 

We also confirmed the criticality of the self-dual point ([8]) of generator II on hypergraph 
A by a Monte-Carlo simulation of 150000 samples, and find results similar to Fig. [8] but 
will not be shown here. In addition, we carried out extensive simulations for the case of 
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Figure 7. Replacing the two central parallel bonds (dashed lines) having probability 
p by one bond having effective probability p' = 2p — p 2 , necessary when Generator I 
on hypergraph A is used. The dotted lines represent the interior of the generator with 
probabilities defined in Fig. UK 
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Figure 8. Plot of s r ~ 2 P> s vs. s a for the Monte-Carlo simulation of generator I on 
hypergraph A, for the first self-dual point, eqn. (|9| (middle), for p + 0.0001 and the 
same r and t as in ([9| (upper), and for p — 0.0001 and again the same r and t (lower). 
The horizontal plot for the first case confirms that the system is at the critical point 
with those probabilities. 


uniform probabilities with this generator and hypergraph, and found the critical value 
p — r — t — 0.441374 ± 0.000001, which falls well within the bounds given by Wierman 






Percolation on hypergraphs with four-edges 


10 


et al. [ 1 5] Jf] 

0.424072 < p c < 0.463661 (14) 

This is a non-self-dual percolation threshold, specific to the uniform generator II on 
hypergraph A. At this value of p, the connection probabilities are P\ = 0.174841, 
P 2 = 0.0645187, P 3 = 0.0786097, P 4 = 0.153594, P 5 = 0.0375026, and P 6 = 0.0120231. 
The duality conditions (|4]-[6| are far from being satisfied by these values, and the values 
of the p are far from those of the other generators listed in Table [!} 

To find this result we simulated up to 10 8 samples for each value of p using a 
smaller cutoff of 2 15 = 32768. Our Monte-Carlo results are plotted in Fig. |9j Here 
we show an alternate way of analyzing the data, where we plot s T ~ 2 P> s vs. s~ n where 
fl = 72/91 [2H [22] represents the finite-size scaling for smaller s. At p c , we have 
s r ~ 2 P> s = A + so at p c our plot should yield a straight line. This form of plotting 

the data is useful when p is very close to p c . 
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Figure 9. Plot of s r ~ 2 P> s vs. s~ n with Q = 72/91 for Monte-Carlo results of the 
uniform probability generator II p = r = t on hyper graph A. The dashed line represents 
the linear behavior which is expected for p = p c , and suggests p c = 0.441374T0.000001. 
Some representative error bars are shown. 
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° 0.441385 
a 0.441375 
0 0.441372 
x 0.441370 
<>0.441365 


5.3. Generator I on hypergraph B 


Here, we have the two additional bonds, one with probability p\ and the other with 
probability p 2 = 1 — pi- However, since we are looking at self-dual square generators, 


X We mention that there appears to be a minor calculation error in ns]. In section 4.2, which 
concerns the model considered here, the authors give the fourth probability as 0.432569051787763. 
However, using Mathematica to solve the fifth equation in their Table 4, P p [C\] + ^[Cy + P p [Cs] = 
Pp[C 4 ] + Pp[Cs] + P p [Cq], we find instead p = 0.4444066898118085. (The fifth probability they give is 
the solution to the fourth equation in Table 4.) However, because this value falls within their other 
values, it does not alter the final bounds (14). 
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which do not depend on pi, we expect that by changing pi nothing changes, and this 
is exactly what is observed. We performed the simulation for four different values of pi 
(0.00, 0.30, 0.50, and 0.85) and for the three self-dual points given in equations @0. 
the system was found to be critical in all of these cases within numerical accuracy. 

5.4. Generator II on hypergraph B 

We verified that the self-dual point ([8]), is the critical point for the four values of p\ 
used above, consistent with our expectations. We also tested the uniform probability 
critical point (p = r = t = 0.441374) that we found for generator II hypergraph A, and 
observed that in this case, as expected, p\ does matter. The system is at criticality only 
for the two points pi ~ 0.7665 and pi ~ 1 — 0.7665 = 0.2335. 

6. Further analysis of critical manifolds 

For generator I, we further studied the self-dual relations. We independently evaluated 
each of the three conditions given in eqns. © and (|6]) and found p in terms of r 
and t. In this way, we found three equations for p(r\ t ), one for each duality condition: 


Pi(p,r,t) = Pi(p, r, t) =*■ p m (r,t) 
Pi(p,r,t) = P 3 (p,r,t) => p <2> (r,t) 
Ps(p,r,t) = P 6 (p,r,t) =t-p (3) (r,t) 


(15) 

(16) 
(17) 


Figures [lOj [Tl] and [12] show the manifolds for the three self-duality conditions. These 
three manifolds intersect at three points—precisely our self-dual points. 


o.o 


0.5 


1.0 


i-° r 


p 



1.0 


Figure 10. p^\r,t) satisfying P^ m F\. The column at the corner of the plot is the 
artifact of the plotting program caused by a singularity. 
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o.o 



Figure 11. p^ 2 \r,t) satisfying P 2 = P3. 


0.0 



Figure 12. p( 3 \r,t) satisfying P~> = P(,. 

To see the intersection of these surfaces, we set pairs of the p( n \r, t) equations equal 
to each other and find t as a function of r. This gives us three non-independent relations 
t(r): 


p ^ (r, t) = p ^ (r, t) => t ^ (r) 

(18) 

p^\r,t) = p^ 3 \r,t ) =>- t^ 2 \r) 

(19) 

p^ 2 \r,t) —p^\r,t) =>- t^ 3 \r) 

(20) 


where t^(r) is satisfied if both self-duality conditions ([ 4 ]) and ([ 5 ]) are satisfied, and so 
forth. The interesting fact is that if we plot the curves of t ®(r) and t^(r), they 

are very nearly identical in a wide region (Fig. [13]), which implies that the points along 
these curves should be quite close to self-duality. Finally, we plot their differences in 
Fig. [L4j which shows clearly that the curves cross at the three self-dual points. The 
curves remain within 0.00001 of each other over a wide range of r, meaning that we 
“almost” have a manifold of thresholds. 
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Figure 13. Plot of the curves of h 2 \r) and t^\r). There are actually three 

different curves here but they are not distinguishable in this figure. 



Figure 14. 

curves cross 


Plot of the difference between pairs of curves of h n \r) of Fig. 
exactly at self-dual points ([£11). 
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The 


Motivated by the results given in Table [TJ we find Pi(p , r, t) along the curves in Fig. 


self-dual points of generator I (where the three curves cross) along with generator III, 
which is a variant of generator I. It can be seen that the points he almost on a straight 
line; we have no explanation for this surprising result. The point (Pi, P 2 ) f° r generator 
II falls on the same approximate line, however, it is not shown on Fig. [15] as it is farther 
away from other points on the plot. Similar nearly linear curves can be obtained for 


13, P 2 vs. Pi is plotted in Fig. M5l Also on the same plot, we show the values for the 
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Pi 



Figure 15. Plot of P 2 (p, r, t) vs. P\ (p, r, t ) for the points t (l> (r) along the curves in Fig. 
fl3l The circles mark the values of the self-dual probabilities. These are indeed at the 
crossing points of the curves. The black square shows the probabilities for generator 
III, which falls nearly on the curve determined by f (3) (r). 

The same analysis was not possible for generator II because the p( n \r, tfs were 
found to be singular at many points, and solutions for t^ n \r) could not be found. 


7. Results for site percolation 


By choosing a generator in which either all or none of the four vertices are connected 
together, we can generate critical site percolation systems. (For the three-hypergraph 
case where the hypergraph is a regular array of triangles, this procedure yields simply 
the system of site percolation on the triangular lattice.) Thus, we assume P 2 = P 3 = 
P 5 = P e = 0, and the self-duality conditions are satisfied if Pi = P\ = 1/2. Applying 
this to hypergraph A, we find site percolation on the triangular lattice once again (Fig. 


16). In Fig. 17 we consider a stretched version of hypergraph A, which is also self-dual, 


and applying the all-or-none generator to this system yields another lattice with a site 
threshold of 1/2, but with non-planar bonds. 

Applying the all-or-none generator to hypergraph B, we find site percolation on 
the union-jack lattice with non-uniform probabilities pi, (1 — pi) and 1/2, as shown 
in Fig. [18} In fact, as shown in that figure, we can generalize this further by having 
a hypergraph with alternating dual hyperedges, where the blue hyperedge is the dual 
to the red hyperedge, and the system is still self-dual. This yields a generalization 
to a union-jack lattice with probabilities p\ and p 2 as shown in Fig. [TS) ). This is an 
interesting system that interpolates continuously between the uniform union-jack lattice 
(pi = p 2 = 1/2) and the covering lattice or line graph for bond percolation on the square 
lattice (pi = 0 or 1). 
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Figure 16. If each hyperedge is replaced by a site in hypergraph A, we get a triangular 
site percolation system. 



Figure 17. Replacing all hyperedges by a site in this self-dual hypergraph created by 
“stretching out” hypergraph A yields another system with site percolation threshold 
oip = 1/2, but with non-planar crossing bonds (on the left). The average coordination 
number of the site lattice is six. 

7.1. Other lattices 

The finding of an inhomogeneous site percolation model motivated us to look for other 
soluble cases of inhomogenous site percolation models to compare with. Here we list 
several of them. For comparison, we can think of our union-jack system as having sites 
with probabilities p\ , p 2 , P 3 , and p 4 , with pi + p 3 = 1 and p 2 + Pa = 1, implying the 
condition (1 — pi — p 3 )(l — p 2 — p,\) = 0 or 

1 - Pi ~ P 2 ~ P3 ~ P4 + P 1 P 2 + PiP4 + P 2 P 3 + P3P4 = 0 (21) 

For inhomogeneous site percolation on the kagome lattice—the covering lattice of 
the honeycomb lattice—Sykes and Essam |Ej showed that the threshold is given by 

1 - P 1 P 2 ~ P 1 P 3 ~ P 2 P 3 + P 1 P 2 P 3 = 0 (22) 

In the limit p 3 —> 1, this gives pi + p 2 = 1, in which case the system becomes the non- 
planar square-lattice covering lattice, similar to what happens to our union-jack lattice 
when pi —> 0 or 1. 

We can find the site threshold for the inhomogeneous martini lattice [I23j by 
considering a martini generator with the three outer bonds having probability pi, p 2 
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(b) 


Figure 18. Hypergraph B and its equivalent site percolation system, (a) We replace 
the bonds with probability p\ and 1 — pi with green and orange sites respectively. 
More generally, we can also have hyperedges with alternate probability P 2 and 1 — p 2 , 
a configuration which we demonstrate by blue and red squares respectively, (b) 
generalized site percolation on the union-jack lattice with site probabilities pi, 1 — p \, 
P 2 , and 1 -p 2 - 


and p 3 , and triangular part correlated with all vertices connected with probability p\ 
and disconnected with probability 1 — p^. The covering lattice to this is the martini site 
lattice, and using ([!]) we find that the threshold is given by 

1 - P 1 P 2 PA ~ PiPm ~ P 2 P 3 P 4 + P 1 P 2 P 3 P 4 = 0 (23) 

For the covering lattice of the inhomogeneous martini lattice, shown in Fig. [I9]b, the 
site threshold is [2H US ESI EE] 

1 - p x p 2 r 3 ~ P 2 Pzri - P\P 3 r 2 ~ PiP 2 nr 2 - PiP 3 r x r 3 - p 2 p 3 r 2 r 3 
+ PiP 2 P 3 rir 2 + Pip 2 p 3 r \r 3 + pip 2 p 3 r 2 r 3 + pip 2 r 1 r 2 r 3 
+ PiP 3 rir 2 r 3 + p 2 p 3 r x r 2 r 3 - 2p x p 2 p 3 rir 2 r 3 = 0 (24) 

which interpolates between site percolation on the inhomogeneous triangular covering 
lattice (by setting all p’s equal to 1) and site percolation on the inhomogeneous kagome 
lattice (setting all r’s equal to 1). 

While there evidently exists several examples of site thresholds for non-uniform 
systems, none of these yield the union-jack lattice in any limit, and consequently do not 
yield our result given above. 

8. Conclusions 

The self-duality and criticality of a hypergraph with four-edges is determined by the 
duality conditions Q, 0 and 0 given by Bollobas and Riordan. While it is easy 
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Pi 




Figure 19. (a) The inhomogeneous martini lattice. The outer star of bonds of the 
generator has probabilities pi, p -2 and p:i , and the inner triangle bonds have probabilities 
ri, r '2 and r 3 . The covering lattice transformation is also shown on the top generator 
where each bond is replaced with a site having the same probability. Two sites are 
connected by a dashed line if the corresponding bonds are connected, (b) the equivalent 
covering lattice, which is also called the (1 x 1 ) : (2 x 2 ) subnet lattice [25]. 


to construct generators with correlated bonds that satisfy these conditions, it is not 
obvious that one can do it with independently occupied bonds, or independent bonds and 
internal sites. In this paper we showed that this can indeed be done, making use of three 
different four-vertex generators (I, II, and III). We imposed the self-duality conditions 
and found the percolation thresholds for these three generators, as given in ([8-12 ). These 
critical probabilities are independent of the hypergraph, as long as the hypergraph is 
self-dual. Then we carried out Monte-Carlo simulations for generators I and II on the 
two self-dual hypergraphs (A and B) and observed these critical probabilities indeed give 
the percolation threshold for the corresponding generator, independent of the irrelevant 
Pi for hypergraph B. Furthermore, we found for generator II on hypergraph A the 
homogeneous threshold p = r = t = 0.441374 ±0.000001, in agreement with the bounds 
(14) found by Wierman et al. lEf Because the uniform generator is not self-dual, this 
homogeneous threshold is specific to that hypergraph (A), and can result in criticality in 
hypergraph B only if we choose p\ ~ 0.7665 or pi ~ 0.2335. We observe that although 
for generator I there are three self-dual critical points, there are also infinitely many 
points very close to being critical. This can be seen from Fig. [l3j [l4|and[T5| 

Using the all-or-none generator, we can find models for site percolation based 
upon the hypergraphs. With this generator on hypergraph B we found a non-uniform 
union-jack model that interpolates between the uniform union-jack lattice and bond 
percolation on the square lattice, providing a range of simple fully-triangulated systems 
that maybe useful for criticality studies. For example, one can see how the lattice metric 
factor 123 or the number of clusters per site [2£|J varies as one goes from one system to 
the other. 

This work shows that while four-hypergraphs are perhaps more restricted than 
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the case of three-hypergraphs (where many self-dual configurations and also uniform- 
probability generators can easily be found), there are still exact systems with 
independent bond and site occupancies that can be constructed. It would be interesting 
to find other examples of self-dual four-hypergraphs and study other generators than 
the ones considered here, perhaps some that satisfy the duality conditions with two or 
one distinct probability. 
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9. Exact connection polynomials 

Following are exact enumeration polynomials for Pi(p,r,t) through P e (p,r,t), defined 
in Eq. (2) of the paper, for generators I, II, and III, where q = 1 — p, s = 1 — r and 
u = 1 — t: 


9 .1. Generator I 

Pi = 4 q 4 rs 3 t 4 + q 4 s 4 t 4 + 16 q 4 rs 3 t 3 u + 4 q 4 s 4 t 3 u + 20 q 4 r 2 s 2 t 2 u 2 + 24 q 4 rs 3 t 2 u 2 + 6 q 4 s 4 t 2 u 2 + 
8q 4 r 3 stu 3 +20q 4 r 2 s 2 tu 3 +16q 4 r s 3 tu 3 +4q 4 s 4 tu 3 +q 4 r 4 u 4 +4q 4 r 3 su 4 +6q 4 r 2 s 2 u 4 +4q 4 r s 3 u 4 + 

4 4 4 

q s u 


P ‘2 = pq 3 r 2 s 2 t 4 + q 4 r 2 s 2 t 4 + 4 pq 3 rs 3 t 4 +pq 3 s 4 t 4 + 4 pq 3 r 2 s 2 t 3 u + 4q 4 r 2 s 2 t 3 u +1 6pq 3 rs 3 t 3 u + 
4pq 3 s 4 t 3 u + 4 pq 3 r 3 st 2 u 2 + 4q 4 r 3 st 2 u 2 + 23 pq 3 r 2 s 2 t 2 u 2 + 3 q 4 r 2 s 2 t 2 u 2 + 24pq 3 rs 3 t 2 u 2 + 
6 pq 3 s 4 t 2 u 2 +pq 3 r 4 tu 3 + q 4 r 4 tu 3 + 10 pq 3 r 3 stu 3 + 2 q 4 r 3 stu 3 + 21 pq 3 r 2 s 2 tu 3 + q 4 r 2 s 2 tu 3 + 
16 pq 3 rs 3 tu 3 + 4pq 3 s 4 tu 3 + pq 3 r 4 u 4 + 4pq 3 r 3 su 4 + 6 pq 3 r 2 s 2 u 4 + 4pq 3 rs 3 u 4 + pq 3 s 4 u 4 


P 3 = p 2 q 2 r 3 st 4 + 2 pq 3 r 3 st 4 + q 4 r 3 st 4 + 3 p 2 q 2 r 2 s 2 t 4 + 4pg 3 r 2 s 2 t 4 + 4 p 2 q 2 rs 3 t 4 +p 2 q 2 s 4 t 4 + 
4 p 2 q 2 r 3 st 3 u + 8 pq 3 r 3 st 3 u + 4 q 4 r 3 st 3 u + 12 p 2 q 2 r 2 s 2 t 3 u + 16 pq 3 r 2 s 2 t 3 u + 1 6p 2 q 2 rs 3 t 3 u + 
4p 2 q |2 s 4 t 3 n+p 2 q |2 r 4 t 2 n 2 +2pg 3 r 4 t 2 n 2 +q ,4 r 4 t 2 n 2 +llp 2 q ,2 r 3 st 2 n 2 +14p(3 ,3 r 3 st 2 n 2 +g 4 r 3 st 2 n 2 + 
28 p 2 q 2 r 2 s 2 t 2 u 2 + 10 pg 3 r 2 s 2 t 2 n 2 + 24 p 2 q 2 rs 3 t 2 u 2 + 6 p 2 q 2 s 4 t 2 u 2 + 2 p 2 q 2 r 4 tu 3 + 2 pq 3 r 4 tu 3 + 
12 p 2 q 2 r 3 stu 3 + 4pq 3 r 3 stu 3 + 22 p 2 q 2 r 2 s 2 tu 3 + 2 pq 3 r 2 s 2 tu 3 + 16 p 2 q 2 rs 3 tu 3 + 4p 2 q 2 s 4 tu 3 + 
p 2 q 2 r 4 u 4 + 4p 2 q 2 r 3 su 4 + 6 p 2 q 2 r 2 s 2 u 4 + 4 p 2 q 2 rs 3 u 4 + p 2 q 2 s 4 u 4 


P 4 = p 4 r 4 t 4 +4p 3 q , r 4 t 4 +6p 2 q ,2 r 4 t 4 +4pq ,3 r 4 t 4 +q ,4 r 4 t 4 +4p 4 r 3 st 4 +16p 3 q , r 3 st 4 +20p 2 q ,2 r 3 st 4 + 
8 pq 3 r 3 st 4 + 6 p 4 r 2 s 2 t 4 + 24p 3 qr 2 s 2 t 4 + 20 p 2 q 2 r 2 s 2 t 4 + 4p 4 rs 3 t 4 + 16p 3 grs 3 t 4 + p 4 s 4 t 4 + 
4p 3 (/s 4 t 4 + 4 p 4 r 4 t 3 u + 16p 3 (7r 4 t 3 n + 24 p 2 q 2 r 4 t 3 u + 16 pq 3 r 4 t 3 u + 4 q 4 r 4 t 3 u + 16 p 4 r 3 st 3 u + 
64 p 3 qr 3 st 3 u + 80 p 2 q 2 r 3 st 3 u + 32 pq 3 r 3 st 3 u + 24 p 4 r 2 s 2 t 3 u + 9 6 p 3 qr 2 s 2 t 3 u + 80 p 2 q 2 r 2 s 2 t 3 u + 
16p 4 rs 3 t 3 n+64p 3 qrs 3 t 3 n+4p 4 s 4 t 3 n+16p 3 q i s 4 t 3 n+6p 4 r 4 t 2 n 2 +24p 3 (3 , r 4 t 2 n 2 +30p 2 (7 2 r 4 t 2 r( 2 + 
12pq 3 r 4 t 2 u 2 + 24 p 4 r 3 st 2 u 2 + 96 p 3 qr 3 st 2 u 2 + 84 p 2 q 2 r 3 st 2 u 2 + 8 pq 3 r 3 st 2 u 2 + 3 6 p 4 r 2 s 2 t 2 u 2 + 
144p 3 gr 2 s 2 t 2 « 2 + 52 p 2 q 2 r 2 s 2 t 2 u 2 + 24 p 4 rs 3 t 2 u 2 + 9 6 p 3 qrs 3 t 2 u 2 + 6 p 4 s 4 t 2 u 2 + 24 p 3 qs 4 t 2 u 2 + 
4p 4 r 4 tu 3 + 16p 3 qr 4 tn 3 + 12 p 2 q 2 r 4 tu 3 + 16 p 4 r 3 stu 3 + 64 p 3 qr 3 stu 3 + 24 p 2 q 2 r 3 stu 3 + 
24 p 4 r 2 s 2 tu 3 + 9 6 p 3 qr 2 s 2 tu 3 + 12 p 2 q 2 r 2 s 2 tu 3 + 16p 4 rs 3 hu 3 + 64 p 3 qrs 3 tu 3 + 4p 4 s 4 hu 3 + 
16 p 3 ( 7 s 4 tn 3 + p 4 r 4 n 4 + 4 p 3 qr 4 u 4 + 4p 4 r 3 su 4 + 16 p 3 qr 3 su 4 + 6 p 4 r 2 s 2 u 4 + 24p 3 qr 2 s 2 u 4 + 
4p 4 rs 3 u 4 + 16p 3 qrs 3 u 4 + p 4 s 4 n 4 + 4p 3 gs 4 n 4 


-P5 = q 4 r 2 s 2 t 4 + 4q 4 r 2 s 2 t 3 u + 2 q 4 r 3 st 2 u 2 + 2q 4 r 2 s 2 t 2 u 2 


Pq = 2p 2 q 2 r 2 s 2 t 4 + 2 pq 3 r 2 s 2 t 4 + 4p 2 q 2 rs 3 t 4 + p 2 q 2 s 4 t 4 + 8 p 2 q 2 r 2 s 2 t 3 u + 8 pq 3 r 2 s 2 t 3 u + 
16 p 2 q 2 rs 3 t 3 u + 4p 2 q 2 s 4 t 3 u + p 2 q 2 r 4 t 2 u 2 + 2pq 3 r 4 t 2 u 2 + q 4 r 4 t 2 u 2 + 8 p 2 q 2 r 3 st 2 u 2 + 
Spq 3 r 3 st 2 u 2 + 26 p 2 q 2 r 2 s 2 t 2 u 2 + 6 pq 3 r 2 s 2 t 2 u 2 + 24p 2 q 2 rs 3 t 2 u 2 + 6 p 2 q 2 s 4 t 2 u 2 + 2p 2 q 2 r 4 tu 3 + 
2 pq 3 r 4 tu 3 + 12 p 2 q 2 r 3 stu 3 + 4 pq 3 r 3 stu 3 + 22 p 2 q 2 r 2 s 2 tu 3 + 2 pq 3 r 2 s 2 tu 3 + 1 6p 2 q 2 rs 3 tu 3 + 
4 p 2 q 2 s 4 tu 3 + p 2 q 2 r 4 u 4 + 4p 2 g 2 r 3 sn 4 + 6 p 2 q 2 r 2 s 2 u 4 + 4p 2 q 2 rs 3 u 4 + p 2 q 2 s 4 u 4 
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9.2. Generator II 

P\ = 4p 2 q G r 4 t 4 + 8 pq 7 r 4 t 4 + q 8 r 4 t 4 + 16 p 2 q 6 r 3 st 4 + 32 pq 7 r 3 st 4 + 4 q 8 r 3 st 4 + 
24 p 2 q 6 r 2 s 2 t 4 + 48 pq 7 r 2 s 2 t 4 + 6 q 8 r 2 s 2 t 4 +16p 2 q 6 rs 3 t 4 + 32 pq 7 rs 3 t 4 + 4 q 8 rs 3 t 4 + 4 p 2 q 6 s 4 t 4 + 
8 pq 7 s 4 t 4 + q 8 s 4 t 4 + 16p 2 g 6 r 4 £ 3 'U + 32 pq 7 r 4 t 3 u + 4 q 8 r 4 t 3 u + 64 p 2 q 6 r 3 st 3 u + 128 pq 7 r 3 st 3 u + 
16q 8 r 3 st 3 u+16p 3 q 5 r 2 s 2 t 3 u+136p 2 q 6 r 2 s 2 t 3 u+192pq 7 r 2 s 2 t 3 u+24q 8 r 2 s 2 t 3 u+32p 3 q 5 rs 3 t 3 u+ 
144 p 2 q 6 rs 3 t 3 u + 128 pq 7 rs 3 t 3 u + 16g 8 rs 3 t 3 w + 16p 3 g 5 s 4 t 3 « + 5 6 p 2 q 6 s 4 t 3 u + 32 pq 7 s 4 t 3 u + 
4g 8 s 4 t 3 « + 24 p 2 q 6 r 4 t 2 u 2 + 48pq 7 r 4 t 2 u 2 + 6 q 8 r 4 t 2 u 2 + 96 p 2 q 6 r 3 st 2 u 2 + 192 pq 7 r 3 st 2 u 2 + 
24 q 8 r 3 st 2 u 2 + 4p 4 q 4 r 2 s 2 t 2 u 2 + 80 p 3 q 5 r 2 s 2 t 2 u 2 + 320 p 2 q 6 r 2 s 2 t 2 u 2 + 28 8pq 7 r 2 s 2 t 2 u 2 + 
36 q 8 r 2 s 2 t 2 u 2 + 20 p 4 q 4 rs 3 t 2 u 2 + 17 6 p 3 q 5 rs 3 t 2 u 2 + 380 p 2 q 6 rs 3 t 2 u 2 + 192pg 7 rs 3 t 2 'U 2 + 
24g 8 rs 3 £ 2 « 2 + 16p 4 g 4 s 4 t 2 « 2 + 112p 3 g 5 s 4 t 2 « 2 + 124p 2 g 6 s 4 t 2 « 2 + 48 pq 7 s 4 t 2 u 2 + 6 q 8 s 4 t 2 u 2 + 
16p 2 g 6 r 4 tw 3 + 32 pq 7 r 4 tu 3 + 4g 8 r 4 7u 3 + 64 p 2 q 6 r 3 stu 3 + 128pq 7 r 3 stu 3 + 1 6q 8 r 3 stu 3 + 
8 p 4 q 4 r 2 s 2 tu 3 + 128 p 3 q 5 r 2 s 2 tu 3 + 3 68 p 2 q 6 r 2 s 2 tu 3 + 192 pq 7 r 2 s 2 tu 3 + 24 q 8 r 2 s 2 tu 3 + 
64p 4 q 4 rs 3 tu 3 +320p 3 q 5 rs 3 tu 3 +336p 2 q 6 rs 3 tu 3 +12 8 pq 7 rs 3 tu 3 + 16q 8 rs 3 tu 3 +64p 4 q 4 s 4 tu 3 + 
128p 3 q ,5 s 4 tM 3 + 9 6 p 2 q 6 s 4 tu 3 + 32 pq 7 s 4 tu 3 + 4 q 8 s 4 tu 3 + 4p 2 g 6 r 4 'U 4 + 8 pq 7 r 4 u 4 + q 8 r 4 u 4 + 
16 p 2 q 6 r 3 su 4 + 32 pq 7 r 3 su 4 + Aq 8 r 3 su 4 + 2 p 4 q 4 r 2 s 2 u 4 + 32p 3 q 5 r 2 s 2 u 4 + 92 p 2 q & r 2 s 2 u 4 + 
48 pq 7 r 2 s 2 u 4 + 6 q 8 r 2 s 2 u 4 + 16p 4 q 4 rs 3 u 4 + 80 p 3 q 5 rs 3 u 4 + 84 p 2 q 6 rs 3 u 4 + 32 pq 7 rs 3 u 4 + 
4 q 8 rs 3 u 4 + 16 p 4 q 4 s 4 u 4 + 32 p 3 q 5 s 4 u 4 + 24p 2 q 6 s 4 u 4 + 8 pq 7 s 4 u 4 + q 8 s 4 u 4 


P ‘2 = p 4 q 4 r 4 t 4 + 4p 3 q 5 r 4 t 4 + 4p 2 g 6 r 4 t 4 + 4p 4 (7 4 r 3 st 4 + 16p 3 g 5 r 3 st 4 + 16 p 2 q 6 r 3 st 4 + 
Qp 4 q 4 r 2 s 2 t 4 + 24p 3 q 5 r 2 s 2 t 4 + 24 p 2 q 6 r 2 s 2 t 4 + 4 p 4 q 4 rs 3 t 4 + 16 p 3 q 5 rs 3 t 4 + 1 6p 2 q 6 rs 3 t 4 + 
p 4 q 4 s 4 t 4 + 4 p 3 q 5 s 4 t 4 + 4p 2 q 6 s 4 t 4 + 4p 4 q ,4 r 4 t 3 'U + 16p 3 q ,5 r 4 t 3 'U + 16p 2 g 6 r 4 t 3 'U +1 6p 4 q 4 r 3 st 3 u + 
64 p 3 q 5 r 3 st 3 u + 64 p 2 q 6 r 3 st 3 u + 2p 5 q 3 r 2 s 2 t 3 u + 3 6 p 4 q 4 r 2 s 2 t 3 u + 110p 3 q 5 r 2 s 2 t 3 u + 
90 p 2 q e r 2 s 2 t 3 u + 4p 5 q 3 rs 3 t 3 u + 40 p 4 q 4 rs 3 t 3 u + 92p 3 q 5 rs 3 t 3 u + 52 p 2 q 6 rs 3 t 3 u + 2 p 5 q 3 s 4 t 3 u + 
lQp 4 q 4 s 4 t 3 u + 30 p 3 q 5 s 4 t 3 u + 10 p 2 q 6 s 4 t 3 u + 6 p 4 q 4 r 4 t 2 u 2 + 24 p 3 q 5 r 4 t 2 u 2 + 24 p 2 q 6 r 4 t 2 u 2 + 
24 p 4 q 4 r 3 st 2 u 2 + 96 p 3 q 5 r 3 st 2 u 2 + 96 p 2 q 6 r 3 st 2 u 2 + 10 p 5 q 3 r 2 s 2 t 2 u 2 + 92 p 4 q 4 r 2 s 2 t 2 u 2 + 
20 6 p 3 q 5 r 2 s 2 t 2 u 2 + 118 p 2 q 6 r 2 s 2 t 2 u 2 + 24 p 5 q 3 rs 3 t 2 u 2 + 130 p 4 q 4 rs 3 t 2 u 2 + 186 p 3 q 5 rs 3 t 2 u 2 + 
54 p 2 q 6 rs 3 t 2 u 2 + 14 p 5 q 3 s 4 t 2 u 2 + 59 p 4 q 4 s 4 t 2 u 2 + 44 p 3 q 5 s 4 t 2 u 2 + 9p 2 q 6 s 4 t 2 u 2 + 4 p 4 q 4 r 4 tu 3 + 
16 p 3 q 5 r 4 tu 3 + 16 p 2 q 6 r 4 tu 3 + 16p 4 q 4 r 3 stu 3 + 64p 3 q 5 r 3 stu 3 + 64p 2 q e r 3 stu 3 + 16 p 5 q 3 r 2 s 2 tu 3 + 
112 p 4 q 4 r 2 s 2 tu 3 + 192 p 3 q 5 r 2 s 2 tu 3 + 5 6 p 2 q 6 r 2 s 2 tu 3 + 48 p 5 q 3 rs 3 tu 3 + 16 8 p 4 q 4 rs 3 tu 3 + 
120p 3 (7 5 rs 3 7u 3 + 24 p 2 q 6 rs 3 tu 3 + 32 p 5 q 3 s 4 tu 3 + 48 p 4 q 4 s 4 tu 3 + 24 p 3 q 5 s 4 tu 3 + 4p 2 q 6 s 4 tu 3 + 
p 4 q 4 r 4 u 4 +4p 3 q 5 r 4 u 4 +4p 2 q 6 r 4 u 4 +4p 4 q 4 r 3 su 4 + 16p 3 q 5 r 3 su 4 + 16p 2 q e r 3 su 4 +4p 5 q 3 r 2 s 2 u 4 + 
28 p 4 q 4 r 2 s 2 u 4 + 48 p 3 q b r 2 s 2 u 4 + 14p 2 g 6 r 2 s 2 « 4 +12 p b q 3 rs 3 u 4 + 42p 4 q 4 rs 3 u 4 + 30 p 3 q 5 rs 3 u 4 + 
6 p 2 q 6 rs 3 u 4 + 8 p 5 q 3 s 4 u 4 + 12 p 4 q 4 s 4 u 4 + 6 p 3 q 5 s 4 u 4 + p 2 q 6 s 4 u 4 


Ps = p Q q 2 r 4 t 4 + 6 p 5 q 3 r 4 t 4 + 12 p 4 q 4 r 4 t 4 + 8 p 3 q 5 r 4 t 4 + 4p 6 q 2 r 3 st 4 + 24 p 5 q 3 r 3 st 4 + 
48 p 4 q 4 r 3 st 4 + 32p 3 q 5 r 3 st 4 + 6 p 6 q 2 r 2 s 2 t 4 + 36 p 5 q 3 r 2 s 2 t 4 + 72 p 4 q 4 r 2 s 2 t 4 + 48 p 3 q 5 r 2 s 2 t 4 + 
4 p 6 q 2 rs 3 t 4 + 24 p 5 q 3 rs 3 t 4 + 48 p 4 q 4 rs 3 t 4 + 32 p 3 q 5 rs 3 t 4 +p 6 q 2 s 4 t 4 + 6 p 5 q 3 s 4 t 4 +12 p 4 q 4 s 4 t 4 + 
8 p 3 q 5 s 4 t 4 + 4p G q 2 r 4 t 3 u + 24 p 5 q 3 r 4 t 3 u + 48 p 4 q 4 r 4 t 3 u + 32 p 3 q 5 r 4 t 3 u + 1 6p 6 q 2 r 3 st 3 u + 
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96 p 5 q 3 r 3 st 3 u + 192 p 4 q 4 r 3 st 3 u + 128 p 3 q 5 r 3 st 3 u + 2 6 p 6 q 2 r 2 s 2 t 3 u + 150 p 5 q 3 r 2 s 2 t 3 u + 
282 p 4 q 4 r 2 s 2 t 3 u + 172 p 3 q 5 r 2 s 2 t 3 u + 20 p 6 q 2 rs 3 t 3 u + 108 p 5 q 3 rs 3 t 3 u + 180 p 4 q 4 rs 3 t 3 u + 
88 p 3 q 5 rs 3 t 3 u + 6 p 6 q 2 s 4 t 3 u + 30 p 5 q 3 s 4 t 3 u + 42p 4 q 4 s 4 t 3 u + 12p 3 q 5 s 4 t 3 u + 6 p 6 q 2 r 4 t 2 u 2 + 
36 p 5 q 3 r 4 t 2 u 2 + 72p 4 q 4 r 4 t 2 u 2 + 48 p 3 q 5 r 4 t 2 u 2 + 24p 6 q 2 r 3 st 2 u 2 + 144 p 5 q 3 r 3 st 2 u 2 + 


28 8p 4 q 4 r 3 st 2 u 2 + 192 p 3 q 5 r 3 st 2 u 2 + 44 p 6 q 2 r 2 s 2 t 2 u 2 + 23 8p 5 q 3 r 2 s 2 t 2 u 2 + 400 p 4 q 4 r 2 s 2 t 2 u 2 + 
200 p 3 q 5 r 2 s 2 t 2 u 2 + 40 p 6 q 2 rs 3 t 2 u 2 + 178 p 5 q 3 rs 3 t 2 u 2 + 218 p 4 q 4 rs 3 t 2 u 2 + 54 p 3 q 5 rs 3 t 2 u 2 + 
I4p 6 q 2 s 4 t 2 u 2 + 4 8p 5 q 3 s 4 t 2 u 2 + 28p 4 q 4 s 4 t 2 u 2 + 4 p 3 q 5 s 4 t 2 u 2 + 4 p 6 q 2 r 4 tu 3 + 24 p 5 q 3 r 4 tu 3 + 
48p 4 (/ 4 r 4 tM 3 +32p 3 g 5 r 4 tw 3 + 16p 6 (7 2 r 3 st'U 3 +96p 5 g 3 r 3 st'U 3 +192p 4 <2 ,4 r 3 st'U 3 +128p 3 <2 ,5 r 3 st'U 3 + 
36 p 6 q 2 r 2 s 2 tu 3 + 176p 5 g 3 r 2 s 2 t« 3 + 236p 4 g 4 r 2 s 2 tM 3 + 5 6p 3 q 5 r 2 s 2 tu 3 + 40 p 6 q 2 rs 3 tu 3 + 
120p 5 g 3 rs 3 7u 3 + 64p 4 g 4 rs 3 £« 3 + 8 p 3 q 5 rs 3 tu 3 + 16p 6 g 2 .s 4 7u 3 + 16p 5 g 3 s 4 7u 3 + 4p 4 g 4 s 4 tw 3 + 


p 6 q 2 r 4 u 4 + 6 p 5 q 3 r 4 u 4 +12 p 4 q 4 r 4 u 4 + 8 p 3 q 5 r 4 u 4 + 4 p 6 q 2 r 3 su 4 + 24 p 5 q 3 r 3 su 4 + 48 p 4 q 4 r 3 su 4 + 
32 p 3 q 5 r 3 su 4 + 9 p 6 q 2 r 2 s 2 u 4 + 44p 5 g 3 r 2 s 2 tf 4 + 59p 4 q 4 r 2 s 2 u 4 + 14 p 3 q 5 r 2 s 2 u 4 + 10p 6 g 2 rs 3 « 4 + 
30 p 5 q 3 rs 3 u 4 + 16 p 4 q 4 rs 3 u 4 + 2p 3 q 5 rs 3 u 4 + 4 p 6 q 2 s 4 u 4 + 4p 5 q 3 s 4 u 4 + p 4 q 4 s 4 u 4 


Pa = p 8 r 4 t 4 + 8 p 7 qr 4 t 4 + 24 p & q 2 r 4 t 4 + 32 p 5 q 3 r 4 t 4 + 16 p 4 q 4 r 4 t 4 + 4 p 8 r 3 st 4 + 32 p 7 qr 3 st 4 + 
96 p 6 q 2 r 3 st 4 + 128 p 5 q 3 r 3 st 4 + 64p 4 q 4 r 3 st 4 + 6 p 8 r 2 s 2 t 4 + 48p 7 <2r 2 s 2 t 4 + 144p 6 g 2 r 2 s 2 t 4 + 
192p 5 g 3 r 2 s 2 t 4 + 96p 4 g 4 r 2 s 2 t 4 + 4p 8 rs 3 t 4 + 32 p 7 qrs 3 t 4 + 96 p 6 q 2 rs 3 t 4 + 12 8p 5 q 3 rs 3 t 4 + 
64p 4 g 4 rs 3 t 4 + p 8 s 4 t 4 + 8 p 7 qs 4 t 4 + 24 p 6 q 2 s 4 t 4 + 32 p 5 q 3 s 4 t 4 + 16p 4 q 4 s 4 t 4 + 4 p 8 r 4 t 3 u + 
32p 7 qr 4 t 3 u + 9 Qp 6 q 2 r 4 t 3 u + 128 p 5 q 3 r 4 t 3 u + Q4p 4 q 4 r 4 t 3 u + 16 p 8 r 3 st 3 u + 12 8p 7 qr 3 st 3 u + 
2>84p & q 2 r 3 st 3 u + 512 p 5 q 3 r 3 st 3 u + 256 p 4 q 4 r 3 st 3 u + 24 p 8 r 2 s 2 t 3 u + 192 p 7 qr 2 s 2 t 3 u + 
564p 6 g 2 r 2 s 2 t 3 « + 720 p 5 q 3 r 2 s 2 t 3 u + 336 p 4 q 4 r 2 s 2 t 3 u + 1 6p 8 rs 3 t 3 u + 128 p 7 qrs 3 t 3 u + 
360 p 6 q 2 rs 3 t 3 u + 41 6p 5 q 3 rs 3 t 3 u + 160 p 4 q 4 rs 3 t 3 u + 4 p 8 s 4 t 3 u + 2>2p 7 qs 4 t 3 u + 84 p 6 q 2 s 4 t 3 u + 
80 p 5 q 3 s 4 t 3 u + 16 p 4 q 4 s 4 t 3 u + 6 p 8 r 4 t 2 u 2 + 48 p 7 qr 4 t 2 u 2 + 144 p 6 q 2 r 4 t 2 u 2 + 192 p 5 q 3 r 4 t 2 u 2 + 
96p 4 q 4 r 4 t 2 u 2 + 24 p 8 r 3 st 2 u 2 + 192 p 7 qr 3 st 2 u 2 + 57 6p 6 q 2 r 3 st 2 u 2 + 76 8p 5 q 3 r 3 st 2 u 2 + 
384p 4 q ,4 r 3 ,st 2 'U 2 + 36 p 8 r 2 s 2 t 2 u 2 + 288 p 7 qr 2 s 2 t 2 u 2 + 812 p G q 2 r 2 s 2 t 2 u 2 + 944 p 5 q 3 r 2 s 2 t 2 u 2 + 
36 8p 4 q 4 r 2 s 2 t 2 u 2 + 24 p 8 rs 3 t 2 u 2 + 192 p 7 qrs 3 t 2 u 2 + 472 p 6 q 2 rs 3 t 2 u 2 + 400p 5 g 3 rs 3 t 2 M 2 + 
64p 4 g 4 rs 3 t 2 « 2 + Qp 8 s 4 t 2 u 2 + 48 p 7 qs 4 t 2 u 2 + 92 p 6 q 2 s 4 t 2 u 2 + 32 p 5 q 3 s 4 t 2 u 2 + 2p 4 q 4 s 4 t 2 u 2 + 
4p 8 r 4 tu 3 + 32p 7 qr 4 tu 3 + 9 6p 6 q 2 r 4 tu 3 + 12 8p 5 q 3 r 4 tu 3 + Q4p 4 q 4 r 4 tu 3 + 16 p 8 r 3 stu 3 + 
12 8p 7 qr 3 stu 3 + 384 p 6 q 2 r 3 stu 3 + 512p 5 q 3 r 3 stu 3 + 25 6p 4 q 4 r 3 stu 3 + 24 p 8 r 2 s 2 tu 3 + 
192p 7 (7r 2 s 2 tM 3 + 496p 6 q ,2 r 2 s 2 t'U 3 + 448p 5 g 3 r 2 s 2 t« 3 + 64p 4 g 4 r 2 s 2 t'U 3 + 16p 8 rs 3 t« 3 + 
12 8p 7 qrs 3 tu 3 + 224 p 6 q 2 rs 3 tu 3 + 64 p 5 q 3 rs 3 tu 3 + 4 p 8 s 4 tu 3 + 32 p 7 qs 4 tu 3 + 1 6p 6 q 2 s 4 tu 3 + 
p 8 r 4 u 4 + 8p 7 qr 4 u 4 + 24 p 6 q 2 r 4 u 4 + 32p 5 q 3 r 4 u 4 + 16 p 4 q 4 r 4 u 4 + 4p 8 r 3 s« 4 + 32 p 7 qr 3 su 4 + 
96 p 6 q 2 r 3 su 4 + 128p 5 g 3 r 3 s« 4 + 64p 4 g 4 r 3 s« 4 + 6p 8 r 2 s 2 M 4 + 48p 7 qr 2 s 2 u 4 + 124p 6 q 2 r 2 s 2 u 4 + 
112p 5 q ,3 r 2 ,s 2 'U 4 + 16p 4 (7 4 r 2 s 2 « 4 + 4p 8 rs 3 « 4 + 32 p 7 qrs 3 u 4 + 56 p 6 q 2 rs 3 u 4 + 16 p 5 q 3 rs 3 u 4 + 
p 8 s 4 u 4 + 8p 7 qs 4 u 4 + 4p 6 q 2 s 4 u 4 


Pr, = p 4 q 4 r 4 t 4 + 4 p 3 q 5 r 4 t 4 + 4p 2 q 6 r 4 t 4 + 4p 4 (7 4 r 3 st 4 + 16p 3 g 5 r 3 st 4 + 16p 2 g 6 r 3 st 4 + 
6p 4 q 4 r 2 s 2 t 4 + 24 p 3 q 5 r 2 s 2 t 4 + 24 p 2 q 6 r 2 s 2 t 4 + 4 p 4 q 4 rs 3 t 4 + 16j9 3 g 5 rs 3 t 4 + 16p 2 g 6 rs 3 t 4 + 
p 4 q 4 s 4 t 4 + 4p 3 q 5 s 4 t 4 + 4p 2 q 6 s 4 t 4 + 4p 4 q 4 r 4 t 3 u + 16 p 3 q 5 r 4 t 3 u + 16p 2 q 6 r 4 t 3 u + 
16 p 4 q 4 r 3 st 3 u + 64 p 3 q 5 r 3 st 3 u + 64p 2 g 6 r 3 s£ 3 « + 2 8p 4 q 4 r 2 s 2 t 3 u + 100p 3 q ,5 r 2 s 2 t 3 M + 
88 p 2 q 6 r 2 s 2 t 3 u + 24 p 4 q 4 rs 3 t 3 u + 72 p 3 q 5 rs 3 t 3 u + 48 p 2 q 6 rs 3 t 3 u + 8p 4 q 4 s 4 t 3 u + 20 p 3 q 5 s 4 t 3 u + 


8 p 2 q b s 4 t 3 u + 6p 4 q 4 r 4 t 2 u 2 + 24p 8 q b r 4 t 2 u 2 + 24 p 2 q b r 4 t 2 u 2 + 24p 4 q 4 r 3 st 2 u 2 + 9 6 p 3 q b r 3 st 2 u 2 + 
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96 p 2 q 6 r 3 st 2 u 2 + 50 p A q A r 2 s 2 t 2 u 2 + 15 6 p 3 q 5 r 2 s 2 t 2 u 2 + 108 p 2 q G r 2 s 2 t 2 u 2 + 46 p A q A rs 3 t 2 u 2 + 
104p 3 q 5 rs 3 t 2 u 2 + 38p 2 q 6 rs 3 t 2 u 2 +16 p A q A s A t 2 u 2 + 16p 3 q 5 s A t 2 u 2 +4p 2 q 6 s A t 2 u 2 +4p A q A r A tu 3 + 
16 p 3 q 5 r A tu 3 + 16 p 2 q 6 r A tu 3 + 16p A q A r 3 stu 3 + Q4p 3 q 5 r 3 stu 3 + 64p 2 q 6 r 3 stu 3 + 44p A q A r 2 s 2 tu 3 + 
112 p 3 q 5 r 2 s 2 tu 3 + 40 p 2 q 6 r 2 s 2 tu 3 + 32p A q A rs 3 tu 3 + 32 p 3 q 5 rs 3 tu 3 + 8 p 2 q 6 rs 3 tu 3 +p A q A r A u A + 
4 p 3 q 5 r A u A + 4p 2 q 6 r A u A + 4 p A q A r 3 su A + 1 6 p 3 q 5 r 3 su A + 1 6 p 2 q 6 r 3 su A + 11 p A q A r 2 s 2 u A + 
28 p 3 q 5 r 2 s 2 u A + 10 p 2 q 6 r 2 s 2 u A + 8 p A q A rs 3 u A + 8 p 3 q 5 rs 3 u A + 2 p 2 q 6 rs 3 u A 


Pq = 2 p 6 q 2 r 2 s 2 t 3 u + 8 p 5 q 3 r 2 s 2 t 3 u + 8 p A q A r 2 s 2 t 3 u + 4 p 6 q 2 rs 3 t 3 u + 1 6p 5 q 3 rs 3 t 3 u + 
16 p A q A rs 3 t 3 u + 2 p 6 q 2 s A t 3 u + 8 p 5 q 3 s A t 3 u + 8p 4 g 4 s 4 f 3 tf + 10 p & q 2 r 2 s 2 t 2 u 2 + 40 p b q 3 r 2 s 2 t 2 u 2 + 
40p 4 g 4 r 2 s 2 f 2 'U 2 + 20 p 6 q 2 rs 3 t 2 u 2 + 68 p 5 q 3 rs 3 t 2 u 2 + 56p 4 g 4 rs 3 f 2 tf 2 + 10 p 6 q 2 s A t 2 u 2 + 


28p b q 3 s A t 2 u 2 + 11 p A q A s A t 2 u 2 + 1 6p b q 2 r 2 s 2 tu 3 + 64 p b q 3 r 2 s 2 tu 3 + 64 p A q A r 2 s 2 tu 3 + 
32p 6 q 2 rs 3 tu 3 + 80p b q 3 rs 3 tu 3 + 32 p A q 4 rs 3 tu 3 + 16 p 6 q 2 s A tu 3 + 16 p b q 3 s A tu 3 + 4 p A q A s A tu 3 + 
4 p 6 q 2 r 2 s 2 u A + 16 p b q 3 r 2 s 2 u A + 16 p A q A r 2 s 2 u A + 8 p e q 2 rs 3 u A + 20 p b q 3 rs 3 u A + 8 p A q A rs 3 u A + 
4 p 6 q 2 s A u A + 4p b q 3 s A u A + p 4 g 4 s 4 n 4 


P.5 1 . Generator III 

For generator III, the probabilities can be found by hand by going through the different 
configurations of the outside bonds first, resulting in the first entry. The second entry 
is the expansion which was verified by exact enumeration: 


Pi = q A (/ (4rs 3 + s 4 ) + u) = 4q A rs 3 t + q A s A t + q A u 


P 2 = pq 3 (t (r 2 s 2 + 4rs 3 + s 4 ) + u)+q A r 2 s 2 t = pq 3 r 2 s 2 t+q A r 2 s 2 t+4pq 3 rs 3 t+pq 3 s A t+pq 3 u 


P 3 = P 2< f (t if s 3 + s) + u) + 2pq 3 rs (1 — s 2 ) t + q A r 3 st = p 2 q 2 r 3 st + 2pq 3 r 3 st + q A r 3 st + 
3 p 2 q 2 r 2 s 2 t + 4 pq 3 r 2 s 2 t + 4 p 2 q 2 rs 3 t + p 2 q 2 s A t + p 2 q 2 u 

P 4 = p A + 4p 3 q + 4 p 2 q 2 r (1 — s 3 )t + 2 p 2 q 2 (1 — s 2 ) 2 1 + 4pq 3 r 2 (1 — s 2 )t + q A r A t = 
p A r A t + 4p 3 qr A t + 6 p 2 q 2 r A t + 4pq 3 r A t + q A r A t + 4p A r 3 st +16 p 3 qr 3 st + 20 p 2 q 2 r 3 st + 8 pq 3 r 3 st + 
6 p A r 2 s 2 t + 24 p 3 qr 2 s 2 t + 20 p 2 q 2 r 2 s 2 t + 4p 4 rs 3 f + 16 p 3 qrs 3 t + p 4 s 4 t + 4 p 3 qs A t + p 4 « + 4p 3 qu 


P5 = q A r 2 s 2 t 


Pq = p 2 q 2 (/ (2r 2 s 2 + 4rs 3 + s 4 ) + «) + 2 pq 3 r 2 s 2 t = 2 p 2 q 2 r 2 s 2 t + 2pq 3 r 2 s 2 t + 4 p 2 q 2 rs 3 t + 
p 2 q 2 s A t + p 2 q 2 u 



